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The parametric deformations of quasienergies and eigenvectors of unitary operators are applied 
to the design of quantum adiabatic algorithms. The conventional, standard adiabatic quantum com- 
putation proceeds along eigenenergies of parameter-dependent Hamiltonians. By contrast, discrete 
adiabatic computation utilizes adiabatic passage along the quasienergies of parameter-dependent 
unitary operators. For example, such computation can be realized by a concatenation of parame- 
terized quantum circuits, with an adiabatic though inevitably discrete change of the parameter. A 
design principle of adiabatic passage along quasienergy is recently proposed: Cheon's quasienergy 
and eigenspace anholonomies on unitary operators is available to realize anholonomic adiabatic al- 
gorithms [Tanaka and Miyamoto, Phys. Rev. Lett. 98, 160407 (2007)], which compose a nontrivial 
family of discrete adiabatic algorithms. It is straightforward to port a standard adiabatic algorithm 
to an anholonomic adiabatic one, except an introduction of a parameter \v), which is available to 
adjust the gaps of the quasienergies to control the running time steps. In Grover's database search 
problem, the costs to prepare \v) for the qualitatively different, i.e., power or exponential, running 
time steps are shown to be qualitatively different. Curiously, in establishing the equivalence between 
the standard quantum computation based on the circuit model and the anholonomic adiabatic quan- 
tum computation model, it is shown that the cost for \v) to enlarge the gaps of the eigenvalue is 
qualitatively negligible. 

PACS numbers: 03.67.Lx, 03.65.Vf 



I. INTRODUCTION 

An adiabatic passage along an eigenenergy of a Hamil- 
tonian with slowly-varying parameter [JJ provides one of 
the simplest ways to control quantum states. Recently, 
one of largest-scale applications of the adiabatic passage 
to compose an algorithm for classically intractable prob- 
lems is proposed by Farhi et al. [2]. Their adiabatic 
passage connects a quantum state, which is supposed 
to be easy to prepare, to a state that represents a so- 
lution of, for example, a satisfiability problem. Farhi 
et al. showed a systematic way to design a parameter- 
dependent Hamiltonian with local interactions of qubits 
only which governs the adiabatic passage. We will call 
this approach the standard adiabatic quantum compu- 
tation (SAQC) as we later introduce a new approach 
to be distinguished from the SAQC. Recently, the stan- 
dard adiabatic quantum computation is proven to have 
the same computational power as the standard quan- 
tum computation in terms of the computational com- 
plexity |3]. It is still an open question whether the adi- 
abatic approach really solve the classically intractable 
problems efficiently as we can see the promising numeri- 
cal results jU [S] as well as the reports of disastrous slow- 
down [sua m. 

In many studies of SAQC, e.g., in the simulations of 
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SAQC by quantum circuits [2] and classical digital com- 
puters 9 , or in experimental realizations where the adia- 
batic changes of coupling constants are infeasible [Till ITT] , 
it is inevitable to introduce the discretization of the adi- 
abatic deformation of parameters. In such cases, the 
time evolution in the computational process is realized 
by the products of parameterized unitary operators, each 
of which represents a computational step to emulate the 
time evolution for the unit of time. Due to the tolerance 
of the adiabatic passage against small perturbations, the 
discretization of the time evolution, in general, provides 
a good approximation of SAQC. The discretization pro- 
cess however can be considered to introduce an alterna- 
tive model of adiabatic quantum computation. The state 
vector follows the adiabatic change of an eigenvector of 
a unitary operator, where the definition of eigenenergy 
can be inapplicable even in any approximate sense. This 
scheme, which we call discrete adiabatic quantum com- 
putation (DAQC), will be formulated based on the adi- 
abatic passage along a quasienergy jT5], or equivalently, 
an eigenangle of slowly- varying unitary operator. 

DAQC is useful to generalize quantum adiabatic al- 
gorithm, which we will show in this paper. A family 
of DAQC that essentially relies on the adiabatic passage 
along a quasienergy is recently proposed by one of the au- 
thors [T5] . Here the adiabatic passage is composed with 
a help of Cheon's eigenvalue anholonomy [T51[T^l[T51ll6j . 
which enables us to design adiabatic passages that visit 
all eigenstates of a given unperturbed Hamiltonian. This 
adiabatic scheme, which will be called anholonomic adi- 
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abatic quantum computation (AAQC), composes an in- 
teresting and nontrivial family of DAQC. In particular, 
AAQC does not approximate SAQC, and hence the ques- 
tion whether AAQC is equivalent to the standard quan- 
tum computation naturally arises |13j . 

In this paper, we establish a formulation of AAQC and 
elucidate its equivalence to the standard quantum com- 
putation. To achieve this, we start from the formulation 
of DAQC, the most general family of the adiabatic al- 
gorithm, in Section [TTJ It turns out that AAQC, which 
composes a subset of DAQC, offers us a system atic way 
to design nontrivial instances of DAQC (Section |lll| . We 
show how the performance of AAQC can be evaluated in 
the Grover's unstructured database search problem [17] . 
It is shown that an ingredient of AAQC strongly af- 
fects the performance that is determined by the "gap" 
of eigenvalues (Section IV I. In Section [V] we show that 
the anholonomic adiabatic quantum computation model 
is equivalent with the standard quantum computation, 
through a modification of Aharonov et al.'s proof [3] of 
the equivalence between the standard quantum compu- 
tation and the standard adiabatic computation model. 
There, it turns out that the control of the gap in AAQC 



discussed in Section IV plays crucial role to show the 
equivalence. 



II. DISCRETE ADIABATIC QUANTUM 
COMPUTATION 



In this section, we will introduce discrete adiabatic 
quantum computation (DAQC) in order to facilitate our 
study of AAQC in the following sections. First, a formu- 
lation of DAQC, which offers a unified framework for the 
adiabatic algorithm, is established in Section |II A| Sec- 
ond, it is shown that a family of DAQC approximates 



SAQC in Section II B Third, we show the equivalence of 
DAQC and the standard quantum computation in Sec- 
tion QTC] 



A. Formulation of DAQC 

We will design a computational system that involves n 
qubits. This naturally introduces the concept of locality 
among qubits in the computational system. Namely, if 
an operation involves only few, n-independent numbers 
of qubits, the operation is called local. The time evolution 
of the quantum state is governed by a unitary operator 
U 8 with an "adiabatic" parameter s. In the stroboscopic 
description of the time evolution, a unit step evolution 
of a quantum state from \tp') to is described by a 
quantum map [18] 



W') = UsW). 



(i) 



be compatible with the one for quantum circuits. First, 
we explain a construction of U s from a product of local 
unitary operations. The unitary operator U s accordingly 
has a natural counterpart of a quantum circuit, where 
each factor of U s corresponds to an element of the circuit. 
If the number of the factors is bounded asymptotically 
by a polynomial of n, U s is called to be efficiently imple- 
mentable. The number of the iterations of the quantum 
map characterizes the computational complexity of the 
algorithm. 

Next we discuss the case that U s is realized by a 
time evolution induced by a time-dependent Hamiltonian 
H s (t) during a finite time interval (say, < t < T), i.e., 



U s = cxp - 



H s (t)dt 



(2) 



Since the quantum map is iterated extensively to realize 
an adiabatic passage, U s is required to be efficiently im- 
plement able, whose definition will be introduced so as to 



where exp<_ denotes the time-orderd exponential. We 
restrict ourselves to the case that H s (t) is composed by 
a finite number of local interaction terms. In particular, 
we call H s (t) efficiently implement able, when the number 
of the local interaction terms for H s (t) is bounded by a 
polynomial of n. Based on the argument by Lloyd [21119]. 
the efficiency discussed here is equivalent to the efficiency 
for quantum circuits mentioned above. In addition to the 
two approaches of the implementation of U s , it is possible 
to merge two approaches to construct U s . An example 
will be shown in Section [V] 

To introduce an adiabatic passage for a quantum map, 
the eigenvalue problem of U s is explained more precisely. 
An eigenvalue of U s is a unimodular complex number, 
due to the unitarity. We parameterize the eigenvalue 
e ~ ie ( s ) by a rea i number 9(s), which is called an eige- 
nangle [20]. By definition, the space of eigenangle has 
a period of 2tt. Since we assumed in the following that 
a unit time interval T is associated with the quantum 
map U S1 e.g. U s is induced by a Hamiltonian time evolu- 
tion during time interval T, an eigenangle 6(s) provides 
a characteristic energy scale E(s) = 8(s)hT~' 1 , which 
is called quasienergy [12] and determined up to mod- 
ulo 2nhT~ 1 . Corresponding to E(s), let |£(s)) be an 
eigenvector of U s , i.e., they satisfy the eigenvalue equa- 
tion U s \£(s)) = e- iTE ^/ H \^(s)). We further require that 
E(s) and |£(s)) are smooth with respect to s. The eigen- 
vector |£(s)) is also called a stationary state, due to its 
characteristic in dynamics. For the fixed value of s, the 
state |£(s)) is invariant under the time evolution induced 
by successive application of U s up to the accumulated 
dynamical phase. 

To explain a discrete adiabatic behavior of the state 
vector, let us change the value of s from to s max , under a 
successive application of U s during L (^> 1) steps. Let si 
be the value of s at l-ih step, where .s = and Sl = s max . 
The exact final state is 

^^^...^(O)). (3) 

The adiabatic theorem for quantum maps [5] ensures that 
converges to |£(s ma x)) up to the global phase as 
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L — > oo, if the gaps between the nearest eigenangles are 
large enough. In reality, L is always kept to be finite and 
the adiabatic time evolution is consequently not exact. 
A criterion to determine L is the condition that the time 
evolution obeys the adiabatic time evolution within an 
error e, i.e., - |£(sm ax ))(£(s max )||| < e. This 

is equivalent with L > L(e), where L(e) is an appropri- 
ate lower bound. Here, the most important source of the 
error in the asymptotic regime L — > oo is nonadiabatic 
transition [21]. Hence the crudest estimate of L(e) will 
be determined by the Landau-Zener-Stiickelberg formula 
with a modification that comes from the use of quasiener- 
gies [25] . instead of eigenenergy, for the discrete adiabatic 
processes. 

To summarize the above argument, a definition of 
discrete adiabatic computation (DAQC) is given: An 
instance of DAQC is specified by an efficiently imple- 
mentable unitary operator U s , an initial state l^o), and 
a schedule {si} l=0 of the adiabatic parameter from sq = 
to Si = s max with the number of steps L. DAQC is de- 
signed to transport the initial state vector |v?o)> which 
is the eigenstate |£(s = 0)) of U , to the target state 
|£(s = s max )) within an error e through the continuous 
deformation of an eigenstate |£(s)) of U s . The initial 
state is be prepared from |0) by an efficiently imple- 
mentable unitary operator £/j, i.e. \^f ) = C7"i |0) . We 
need to make L large enough so that the final state 
\^l) = U SL U SL _ 1 ...U 8l \^o) of the computation is a 
good approximation of the target state of the compu- 
tation |f(Smax))- The running time L{e) is the minimum 
value of L to achieve that the distance between the final 
state I^l) and |£(s max )) is smaller than a given error e. 



part) is 



< M> 1 -i J H(s(t))dt\ |* ) 



(4) 



verges to the target for the final state. The running time 
of the computation is the minimum value of T max that 
keeps the distance between the finial and the target states 
within a given error e. 

A discretization of the time-dependence in s(t) gives 
us an instance of DAQC. For a discretization of time t\ 
(0 < I < L), we introduce a piecewise constant function 
S(t) = ai = [a(tj_i) + s(ti)] /2 for i t _j < t < U. For 
the sake of simplicity, we set sq = and t_\ = —t\. 
Accordingly, we have a unitary operator of a time interval 



Ui = exp -iH(si)(ti-ti-i) 



The final state of a SAQC with the schedule s(t) is 



U L U L . 



(5) 



(6) 



Namely, we obtain a DAQC from a SAQC, since |\&o) 
is also an eigenvector of Uq — e^ lH ^°^ t °~ t - 1 \ We may 
expect the convergence of a DAQC to the original SAQC 
in the limit L — > oo. For finite L, the DAQC is an ap- 
proximation of the SAQC. Hence we may regard that the 
DAQC approximately follows the adiabatic passage built 
on the parametrically varying eigenenergy. 



B. DAQC from SAQC through the discretization 
of continuous-time evolution 

As a way to construct DAQC, the discretizations of 
standard adiabatic quantum computation (SAQC) are 
explained for a pedagogic purpose. As mentioned in the 
introduction, many studies of SAQC are done through 
the discretization of the variation of the adiabatic pa- 
rameter. The discretized time evolution is introduced 
so that it approximates well the ideal, continuous-time 
evolution. We first explain the standard adiabatic quan- 
tum computation. A computational system is specified 
by a Hamiltonian H(s), which can be decomposed into 
a finite number of local Hamiltonians, with an adiabatic 
parameter s. We assume, without loss of generality, that 
the initial and the final values of s are and s max , re- 
spectively. The initial state |^o) must be the ground 
state of Hb = H(0), which can be efficiently prepared. 
The target of the computation is the ground state of 
Hp = H (s max ). For a given time-dependence of s(t) with 
< t < T max , s(0) = 0, and s(T max ) = s max , the exact 
final state of the computation (of the unitary evolution 



C. DAQC is equivalent to the standard quantum 
computation 

Now we examine the equivalence of DAQC and the 
standard quantum computation, in the sense of compu- 
tational complexity. First, it is straightforward to see 
that a DAQC is efficiently simulated by a quantum cir- 
cuit from the definition above. Second, for the inverse 
argument, i.e., whether a quantum circuit can be simu- 
lated by a DAQC, we need to remind that DAQC involves 
a rather wide class of computation. In particular, as is 
explained above, DAQC contains the discretizations of 
the standard adiabatic quantum computation, which are 
proven to be equivalent with the standard quantum com- 
putation [3]. Hence, it is trivial to conclude that DAQC 
are equivalent with the standard quantum computation. 

Much more serious question on the equivalence arises 
when we examine a subclass of DAQC, by introducing 
an alternative design of DAQC, where the adiabatic pas- 
sage built on the parametrically varying eigenenergy is 
inapplicable. Instead, an adiabatic passage built on a 
quasienergy needs to be employed. We will discuss such 
an example in the next section. 
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III. ANHOLONOMIC ADIABATIC QUANTUM 
COMPUTATION 



In this section, anholonomic adiabatic quantum com- 
putation (AAQC), which compose a subclass of DAQC, 
is introduced. The adiabatic passages of AAQC can 
be constructed with the help of Cheon's eigenvalue and 
eigenspace anholonomies, which will be explained in Sec- 



tion III A Subsequently, a formulation of AAQC is intro- 



duced in Section |III B| Finally, the simulation of AAQC 
by quantum circuits is examined in Section [ill C| where a 
crucial ingredient to determine the cost of the simulation 
is identified. 



A. Cheon's anholonomies for unitary operators 

First of all, we explain Cheon's eigenvalue and 
eigenspace anholonomies |15j and their realization in 
a periodically pulsed system under a rank-1 perturba- 
tion [13] • This offers a systematic design principle for adi- 
abatic passages, in particular, adiabatic quantum com- 
putation, along parametric changes of unitary opera- 
tors [16]. We assume that the spectrum of an "un- 
perturbed" Hamiltonian H a is discrete, finite and non- 
degenerate. We will examine a periodically-pulsed driven 
system described by the following "kicked" Hamiltonian: 



H s (t) = H + s\v){v\ S(t - nT), 



(7) 



where T and s are the period and the strength of the 
perturbation, and \v) is assumed to be normalized. The 
kicked Hamiltonian ^ naturally induces a quantum map 
IVVi+i) = U s \t/j n ), where \ip n ) is the state vector just be- 
fore the kick at t = nT and 



U, 



e -iH T e -is\v)(v\ 



(8) 



is the corresponding Floquet operator. We set ft = 1 in 
the following. 

We remark the topological structure of the parameter 
space s for the Floquet operator U s |8]). From e~ ls \ v ^ v \ = 
(1 — + e~ ls \v)(v\ , both U s and its spectrum is 

periodic in s with period 2tt. Hence we may identify the 
parameter space of s as a circle. 

We explain what will happen for each quasienergy and 
each eigenvector of U s |8| when we increase s for a period. 
Let us consider the ground energy E' and the first excited 
eigenenergy E" of Hq (E' < E") [23] • The corresponding 
normalized eigenvectors are denoted by \E') and \E"), re- 
spectively. For simplicity, we employ a normalized vector 
\v) = a\E')+b\E"), where a ^ 0, b ^ and |a| 2 +|6| 2 = 1, 
so that the subspace spanned by \E') and \E") is invari- 
ant under U s . Accordingly, in the subspace, which we 
will focus on, U s is a 2-dimensional matrix effectively 
Let |£o(s)) be the eigenvector of U s with |£o(0)) = \E') 
and Eq(s) the corresponding quasienergy, whose branch 



is chosen as [E',E' +2tvT- 1 ), i.e., E' < E" < E+2-kT- 1 . 
Since d s E Q (s) = T' 1 {£ (s)\v) {v\£ (s)) is positive, E (s) 
increases as s increases. At the same time, because of the 
periodicity of C/ 27r = U , E Q (2ir) = E' + [d s E (s)] ds 
must be equal to E' or E" (modulus 2-itT- 1 ). The 
above choice of \v) ensures < d s Eo(s) < 2irT~ 1 . 
i.e., E' < E (2tt) < E' + 2-irT- 1 . Hence we have 
Eq(2tt) = E" , which implies a presence of quasienergy 
anholonomy. The corresponding eigenvector |£o(s)) ar- 
rives at \E") at s = 2tt. Thus the minimum example 
of a path along quasienergy anholonomy Eq(s) is shown. 
Along the path, an adiabatic increment of s for a period 
of its parameter space (i.e., from to 2n) transfers the 
state vector that is initially prepared to be \E'), to \E"). 

So far, we did not specify \v), except that it is a su- 
perposed state of \E') and \E"), because the quasienergy 
anholonomy do not depend on the details of \v). For 
example, the quasienergy anholonomy persists even un- 
der adiabatically slow fluctuations of \v). On the other 
hand, \v) affects on the precise shape of Eq(s) and ad- 
jacent quasienergy curves. Namely, \v) determines the 
quasienergy gaps between neighboring levels of Eq(s) and 
accordingly a timescale for adiabatic passages. 

Extensions of the quasienergy anholonomy for multi- 
levels are straightforward [13j EE]. For example, if we 
choose \v) that has non-zero overlapping with all eigen- 
vectors of H , E (2Trn) is the n-th excited energy of 
Hq [23]. Thus Eq(s) connects all the eigenvalues of Hq 
to offer a path for adiabatic passages with adiabatic in- 
crements of s. 



B. Formulation of AAQC 

An anholonomic adiabatic quantum algorithm can 
be described based on the standard adiabatic algo- 
rithm. The standard adiabatic algorithm is specified by 
a parameter-dependent Hamiltonian H (s) on the Hilbert 
space Ha- Here, we employ the simplest case that H(s) 
is linear in s: H(s) = (1 — s)H-q + sHp, where Hb and 
Hp are the initial and the final Hamiltonians, respec- 
tively. We impose that Hb has non-degenerate ground 
state |0b) with the ground energy 0. The finial Hamil- 
tonian Hp is a cost function of the problem to solve. 
Namely, an eigenvalue of Hp indicates a "distance" be- 
tween the corresponding eigenvector and the answer of 
the problem. We assume that the answer is unique and 
let \x) be the ground state of Hp. To reuse Hb and Hp 
in the anholonomic approach, the maximum eigenvalues 
of Hb and Hp need to be finite. 

We now explain an implementation of an anholonomic 
adiabatic quantum processor under the assumption that 
the standard adiabatic quantum processor is available. 
The state space Ha, and its Hamiltonians Hb and Hp 
are reused. An additional qubit is employed as a "con- 
trol register," whose Hilbert space He has a complete 
orthonormal system {|0), |1)}. Hence the whole Hilbert 
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space is Ha ® He- Corresponding to the "initial" and 
"final" states of computation, we introduce two projec- 
tors Z^ ee |0)(0| and 1^ ee |1)(1|, respectively, on He- 
Combining these parts, we obtain an unperturbed Hamil- 
tonian 



H = H B ® + (E P + Hp) <g> J (c \ 



(9) 



where we assume that Ep is positive and smaller than 
the first excited energy of Hp, . It is straightforward to see 
that the ground energy and the first excited energy Ep 
of Hq are nondegenerate. By construction, the ground 
and the first excited states of Hq are | — ) = |0b) <8> 
|0) and | + ) ee \x) ® respectively. We employ the 
quantum map (|8| to realize the quasienergy anholonomy 
that connects between the eround energy and the first 
excited energy Ep of Hq (|9|. We impose the conditions 
( ± \v) for \v) so that the quasienergy E (s) of U s 
with Eq(0) — reaches Ep at s — 2tt. The adiabatic 
passage along Eq(s) (0 < s < 2tt) transfers the state 
vector from | — ) to | + ) , then the final state provides the 
answer. To evaluate the performance of the anholonomic 
adiabatic quantum computation, we need to obtain the 
quasienergy gap with respect to the path E (s). The 
quasienergy gap depends on the choice of \v). We will 
evaluate the gaps for two choices of \v) in the following 
sections. 



C. Costs for simulations of AAQC by quantum 
circuits 

We show how the unitary operator U s ^ of the an- 
holonomic adiabatic quantum processor is composed by 
local quantum circuits, in a similar way explained in 
Farhi et al. [2]. From the definition of Hq Q, U s is 
expressed by a product U s = Uqb Uqp e~ IS '^^, where 

Uq B EE e -*H B ®Z^T ^ (j Qp = e -i(ff P +Ep)®£<<=>T 

In the analysis of Uqb, we assume that Hb is composed 
by Hb.j , initial Hamiltonians for j-th qubit. For example, 
it is often the case that Hbj is an Hadamard operation. 
This allows further decomposition: 



iZ {c ''T 



Uqb = e - iH **® z ' a)T e- iH ^ 

(10) 

where e ^H B , 3 ®z^T = (]_ _ jj(c)) + Z^e- iH ^i T is a 
controlled- 1 bit unitary operation. Hence Uqb is simu- 
lated by n controlled-1 bit unitary gates. 

To examine Uqp, assume that Hp is a cost Hamilto- 
nian of a satisfiability problem, which is composed by 
clauses Cj. Let m be the number of the clauses and 
Hp j a corresponding local Hamiltonian for Cj. Note 
that Hp = X^-ffpj, and the elements Hpj commute 
with each other. Hence we have a decomposition of E/qp : 



Uqp = e - l£ p* (c)T e -^p,o®* <c) T 



where e-^^x^T = ^ _f(c)) + j(c) e -tH P , 3 ig a gev _ 
eral qubits unitary operation with a condition (3 qubit 
unitary operation for 3-SAT problem). Hence, a simula- 
tion of Uqp requires a conditional phase-shift gate and m 
conditional, few bits unitary gates. 

Finally, e~ ls ^^ is examined. With two unit vectors 
l^o) and in 7^a, \v) is written as \v) = col^o) ® |0) + 
ci | Vi)<8> 1 1) , where Co and c\ are coefficients. We introduce 
two unitary operators Go and G\ so as to induce \vq) 
and \vx) from |0) S Ha, i-e., \v ) = Gq\0) and \vi) — 
Gt\0). Hence we have \v) = G (|0) ® |C)), where G = 

(G (g)i( c ) +I (c) ) (i< c) + Gi(g)X( c )) and|C) eec o |0) + 
Ci 1 1) , and 



-is\v)(v\ = Q e -ts\0){0\®\C){C\Qj 



(12) 



where e -«|0>(0|®|O><C| ig a p h aS e-shift gate. Note that 
G is a product of two conditional-unitary operations. If 
both Go and G\ are simulated efficiently by quantum 
circuits, in other words, both \vq) and \v\) can be effi- 
ciently prepared by the standard quantum computers, G 
can also be simulated efficiently. For example, a state 
2-«/2 Y^j^Q 1 \ j) in Ha is efficiently available, since the 
state is H® n \0) and is available through n Hadamard 
operations H. We may employ if®"|0) for \vq), and, 
we also employ H® r ' 



in Section 



IV B 



for \vi) . On 

the other hand, in Section IV A| we employ \v\) = \x), 
which is the answer of the problem. Note that if \x) 
is efficiently available with quantum circuits, G is also 
efficiently available with quantum circuits and we may 
put G in our anholonomic processor only with a quali- 
tatively negligible cost. Otherwise, the implementation 
of G requires exponentially many quantum gates by the 
definition of the efficiency here. 



IV. TWO EXAMPLES OF AAQC 

Through the study of two examples, we show that the 
magnitude of the minimum quasienergy gap drastically 
depends on \v). At the same time, we examine the cost 
for the increment of the gap, in terms of the number of 
quantum gates to realize \v). 



A. An "optimal" choice of \v) to widen the 
quasienergy gap 

Let us assume that the unique solution of a given prob- 
lem \x) is available to construct \v) for the anholonomic 
processor: 



1 



V2 



(|Q B ) ® |Q) + \x) <g> |1)) 



V2 



(13) 



x e -i-ff P ,i«>I lc 'T , , , e -iir P , m _iig>Z^T 



(11) 



As is seen below, the point of the assumption is that \v) 
is a superposed state of | ± ) , and the precise values of 
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the coefficients are irrelevant. The following argument is 
rather general since it is applicable as long as the problem 
has a unique solution. 

First, we examine the quasienergy gap during the adi- 
abatic passage. With our choice of \v) (13|, the linear 
span of | ± ) is invariant under U s (Eq. ( 81) , since | ± ) 
are eigenstates of the unperturbed part Uq and the linear 
span of | ± ) is invariant under the kick part e~ ts ^^ of 
due to the fact 



-is\v) {v\ 



1 



1 



(l-ia y ) + e ls -(l + a x ), 



(14) 



Ei(0) 









1+) 


l+> 









Eo(2n) 



where a x = |— )(+|+h.c. and a y = — |+h.c. Hence 
it is suffice to examine U s within the subspace spanned 
by | ± ). The characteristic equation of U s within the 
subspace is 



det 
= 





-iEpT 



-is/2 



COS | 

-i sin ; 



cos 



(15) 



The above characteristic equation is independent of the 
details of Hp, in particular, the size of the problem n. 
This indicates that the quasienergy gap is also indepen- 
dent of n. Hence, the choice of \v) above (13 1 provides 



an extraordinary speed up of the anholonomic processor 
for asymptotically large n. 
The eigenvalues are 



z± (s) = exp 



±6 P (s) 



(16) 



cos cos I 



The correspond- 
s)/2 ± e P (s) are de- 



( 13 1 has to be com- 



where 0p(s) = cos 
ing quasienergies E±(s) = (EpT ■ 
picted in Fig. [T] 

The extraordinary speed up by \v 
pensated by the cost for the preparation of the kick part 
of the Floquet operator e~ ls ^^, as is mentioned in the 
previous section. This means that we have to pay the 
cost to prepare the unitary operator that makes |x) from 
a simple state, say, |0)® n . Namely, the effort to solve the 
given problem by the standard quantum computation is 
required. In this sense, the present choice of \v) is not 
practically useful. However, the present argument is ap- 
plicable to study the theoretical nature of AAQC, as is 
shown in Sec W\ 



B. A "fair" choice of \v) for Grover's problem 

We first explain the standard adiabatic algorithm [5] 
for Grover's quantum search |17j before we examine a 
fair choice of the state |u). Let N be the number of 
items in the unstructured database. The items will be 
labeled with integers 0, . . . , N — 1, and the label of the 
answer is denoted by x, which is assumed to be unique. 



FIG. 1: (Color online) Two quasienergies E±(s) (bold curves) 
for the optimal choice of \v) given at Eq. (| 1 3f > - with T = 1 and 
Ep — 2n/3. The dashed line corresponds to a unperturbed 
eigenenergy -Ei(O) = Ep. The period of quasienergy space is 
Smax = 2-7T. The "ground" quasienergy for the present choice 
of the branch exhibits anholonomy, connecting Eq(s = 0) = 
and Eo(s — 2ir) — Ep. Accordingly the adiabatic passage 
along the ground quasienergy transports | — ) at s — to 
| + ) at s = 2ir. Note that the "first excited" quasienergy that 
connects (s,E) = (0,Ep) and (27r,i? max ) has crossings with 
other quasienergies (not depicted), which are s-independent 
and, whose eigenspaces are orthogonal to | ± ) [TS] . 



The Hilbert space of the arithmetic register for the stan- 
dard adiabatic computation is spanned by an orthonor- 
mal system The finial Hamiltonian Hp is a 

cost function of Grover's problem: 



Hp 



a(l-\x)(x\), 



(17) 



where we introduce an energy scale a(> 0). With an 
initialization of the arithmetic register to be in |0b), an 
adiabatic passage along the path of the ground energy of 
H(s), with an adiabatic change of s, transposes the state 
of the arithmetic register to \x). In SAQC, the mini- 
mum gap for the ground energy of H(s) is 0(N~z) for 
asymptotically large N [2]. This implies that the run- 
ning time is O(N), if ds/dt is kept constant along the 
path. Namely, no quantum speedup is available. How- 
ever, a time-dependent ds/dt introduced by Roland and 
Cerf provides a Grover-type improvement of running time 
0(y/N), once the locations of narrow gaps in the param- 
eter space could be identified [2S]. In the latter sense, 
the standard adiabatic quantum search has the same 
computational power as the original Grover's algorithm. 
We will show that the anholonomic adiabatic approach 
has the same performance as the standard adiabatic ap- 
proach. 

In the Grover's problem, we have no prior knowledge 
of \x) to choose \v). Hence, we employ \T), which need to 
have non-zero overlap with \x), to introduce a normalized 
|«>: 

\v) = o|0b) ® |0) + b\T) <g) |1) = a\ - ) + b\T) ® |1), (18) 

where we choose the phase of \T) so that b is positive. 
Although the overlap (x\J-) is required to be non-zero, 
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KxlJ 7 )! become small as N — > oo, due to the lack of the 
prior knowledge of x to prepare A typical state 

for \T) may be a superposition of all candidates of the 
answer, that is 



JV-l 
VJV TO =0 



(19) 



in this sense we consider the state \T) as a fair choice 
for The preparation of \T) and hence |u) is trivial, 
and hence only polynomially many quantum gates are 
required to prepare U s . Hence it is natural to assume 
(x\F) = 0{N^ 1 / 2 ) as N -> oo. We introduce a small 
parameter e = ©(TV -1 / 2 ) to parameterize the overlap 
between |x) and \T) as (a^J 7 ) = e ie sine, where < e < 
7r/2 and 6 = axg(x\J r ). 

To study the spectrum of U s , it is convenient to intro- 
duce a normalized vector Ix- 1 } £ Ha- 



(l-\x)(x\)\7). 



(20) 



Hence we have an expansion of \v) (18) by orthonormal 
vectors: 

\v) =a\-) + be w sme\+)+bcose\f), (21) 
where a normalized vector 



I/) 



II) eH A ®H c 



(22) 



satisfies ( ± |/) = 0. 

The subspace spanned by |±) and |/) is invariant under 
U s = e -iH T e -is\ v )(v\ (|J of Grover's problem. This is 
because \x ) is an eigenvector of the cost Hamiltonian 
for Grover's problem, Hp (17) and |/) is accordingly an 
eigenvector of H 

Hence it is suffice to examine a truncation of Uq by a 
three-level system 



U = \-)(-\ + e 



-iE P T\ 



•)(■ 



l/X/l (23) 



for Grover's problem. The minimum quasienergy gap 
of U s along an increment of s is immediately obtained 
through numerical diagonalizations of 3 x 3 unitary ma- 
trices (Fig. [2). 

We show that the adiabatic passage from | — ) to | + ) 
encounters a narrow quasienergy gap in Fig. [2] leaving 
the details of calculation to Appendix [X] The magnitude 
of the gap is 0(e) = 0(N~ 1 ^ 2 ) by a perturbation expan- 
sion by e. An application of Landau-Zencr-Stiickelberg 
formula for an estimation of nonadiabatic error implies 
that the running time of adiabatic approach is O(N). 
However, if the location of the quasienergy gap is iden- 
tified, Roland and Cerf's prescription [35] is also appli- 
cable to the quantum map U s to obtain a Grover-type 
improvement of running time 0(y/N). 



E 2 (0) 



Ei(0) 









If > 


If > 




1+) ^ 


l+> 





£i(2tt) 



Eo(2n) 



In 



FIG. 2: (Color online) Three quasienergies E±(s) and Efts) 
(bold curves) for the fair choice of \v) (Eqs. (18 1 and (19l), 
with |(a;|.F>| = 1/100 (i.e., N = 100) and |<"~>>| 2 = 5/6. 



The energy scale for the cost Hamiltonian (17 1 is a = 2tt/3. 
Other parameters are the same as in Fig [l] The adiabatic 
passage for the computation follows the ground quasienergy 
Eo(s), which encounters a narrow avoided crossing during the 
increment of s from to 2n. 



V. EFFICIENT SIMULATION OF QUANTUM 
CIRCUITS BY AAQC 

In this section, we will complete to show the equiva- 
lence of AAQC with quantum circuits. Since we have 
already shown that AAQC is efficiently simulated by 
quantum circuits, we only need to show the converse; 
quantum circuits can be efficiently simulated by AAQC. 
To achieve this, we utilize Aharonov et al.'s SAQC [5] 
to construct a simulator of a given quantum circuit by 
AAQC. We remark that their proof for the equivalence 
of SAQC to quantum circuits can be directory applied to 
DAQC, since a subset of DAQC consists of all discrete 
approximant of SAQC, i.e., SAQC is a subset of DAQC 
in effect. On the other hand, as for AAQC, the relation- 
ship with SAQC is not trivial. This is the reason why the 
equivalence of AAQC to quantum circuits remains non- 
trivial. The quasienergy gap of our simulator is adjusted 
by \v), as is done in Section IV A with a reasonable cost 
for the present case. 



Because the following argument involves extensive use 
of local observables, it is convenient to denote them by 
local operators, instead of brackets. Let Hb denote the 
Hubert space of a qubit and {|0), |1)} corresponding unit 
vectors. Operators of the qubit are 

X= |l)(0| + h.c, Y = i|l)(0| + h.c, 

Z=|0)(0|-|1)<1|, 



i=i(l + i) = |Q)(0|, X: 



1 



(1-Z) = |1}(1|. (24) 



and we need to introduce "annihilation" and "creation" 
operators: 

i=|0)<l|, At = |1><0|. (25) 

The target of the simulation is a quantum circuit that 
involves n qubits, whose Hilbert space is denoted by 
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Ha = {Hh)® n - Let be a local operator a on j-th 
qubit. The quantum circuit is composed by L elemen- 
tary circuits > where ?// involves only few (possi- 
bly two) qubits, and its state evolves as 



\a(l)) = & l \a(l-l)). 



(26) 



whose ground subspace is Ha® |c(0)) , if the state space is 



restricted within Ha < 
circuit, we introduce 



(V) 
clock ' 



To initialize the quantum 



H 



input 



(A) 



2l, 



(31) 



where the initial and final states of the quantum cir- 
cuit are |a(0)) = |0)® n and \a) = \a(L)), respectively. 
Our simulator of the quantum circuit by an AAQC de- 
pends on the initial and final Hamiltonians Hb and Hp 
of Aharonov et al.'s simulator of the same circuit by a 
SAQC. Hence the contsruction of these Hamiltonians ex- 
plained. 

Another ingredient of the simulator is Kitaev's 
clock [25], which is composed by L qubits. Let H c iock 
be the Hilbert space of the clock's states and di be a lo- 
cal operator a on clock's l-th qubit. The state of clock at 
step I is associated with a state 



=(0> 



EE \l l L - 1 ) 



(27) 



which is an eigenstate of X1X2 ■ ■ - T-i^i+i ■ ■ ■ Zl-iZl with 
a non-degenerate eigenvalue 1. Although {\c(l))}f =0 is 
orthonormal, it is not a complete system of H c iock- Let 
us introduce a Hilbert space H^J c y. that is spanned by 
{\c(l)}}f =0 . When the state of the clock does not belong 
to HG2. k , the clock is "out of order". A cost Hamiltonian 



on H 



clock 



L-l 



H, 



clock 



E^* 



1+1 



1=1 



is introduced in order to characterize H 



(V) 
clock' 



(28) 



By dcfini- 



/(V) 



tion, the ground subspace of -ff clock is ?^ lock - 

Combining Kitaev's clock and the quantum circuit, the 
state of the whole system "atZ-thstep" is |7(Z)) = \a(l))(g) 
\c(l)). Let \rf) be a superposition of | / y(i))'s: 



1 



1=0 



(29) 



which is the destination state, i.e., the ground state of 
Hp, of Aharonov et aZ.'s simulator of the quantum cir- 
cuits. Let Ho be L-dimensional Hilbert space spanned 

by {|t(0)}^=o- ^ wm ^ e c l ear that we need to focus on 
Hq in the following analysis of the simulator. 

The initial state of Aharonov et al.'s simulator is 
| 7 (0)) = |0)®(™+ L ). To achieve this, H B , the initial 
Hamiltonian of Aharonov et al.'s SAQC, is composed by 
three parts. The first ingredient is i/ciock (28 1, whose 
ground state must be in Ha <8 H, 



(V) 



clock' 



The next is 



H 



clockinit 



(30) 



whose ground subspace contains (7(0)) and Ha 



(span{|c(Z))} i= 
are combined: 

Hn 



To obtain Hb, three Hamiltonians 



H 



clockinit 



H 



input 



H 



lock; 



(32) 



whose unique ground state is |7(0)). 

The finial Hamiltonian Hp of Aharonov et al.'s SAQC 
is also composed by three parts. Two of them, -ffi npu t and 
H dock, are also contained in Hb- The last part, which 
we denote Hh, of Hp is determined so that \rj) (Eq. (29 1) 
is a nondegenerate ground state of Hp. A component of 
Hh, which correspond to the "Z-th step" of the quantum 
circuit, is a "five-body" interaction term (1 < I < L): 



1 



1 - \ 



_iC/ ;J 4;z i+1 + h. c 



where the "three-body" term 

(UiA 



1 r 

2 



1 



h.c. 



(33) 
(34) 

(35) 



transfers into a superposition of ± 1)) and the 

projectors Ti—i and Zi + \ make Hi to be invariant in Ho- 
With the boundary terms 



H 1 = hxZ 2 , and H L =lL-ihL, 



we have 



Hh = E Hi ; 



(36) 



(37) 



/=i 



which has a unique ground state |ry) with the ground 
energy and leaves Ho invariant. Hence we obtain the 
final Hamiltonian 



Hp 



H h +H 



input 



+ #clock- 



(38) 



In the following, we employ the restriction of Hp to Hq. 

To incorporate Aharanov et al.'s Hp with the AAQC 
simulator, we need to estimate the asymptotic behaviors 
of its energy gap of the ground state and the maximum 
eigenenergy in the limit L —> 00, within the subspace 
Ho- The spectrum properties of Hp are also crucial for 
the Aharonov et al.'s SAQC and are already clarified by 
them [3]. Hence we explain their results only briefly. 
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Since the subspace Ho is an eigenspace of i?i np ut + #ciock, 



which is a part of Hp (see, Eq. ( 38 1 ) , and the correspond- 
ing eigenvalue is zero, it is suffice to examine H^, the 
nontrivial remainder of Hp. It is straightforward to es- 
timate the spectrum of H^ from the fact that H^ is a 
discretized one-dimensional Laplacian. The first excited 
energy of Hp is 



A = <D(L~ 2 ) 



(39) 



and is doubly degenerated. The maximum eigenenergy 
of Hp is 



Wp 



2 + e>(ir 1 ). 



(40) 



Hence the two Hamiltonians Hp, and Hp for Aharonov 
et al.'s SAQC have introduced. 

With these preparations, we introduce an AAQC, 
which simulates the quantum circuit. The construction 
of the AAQC from Aharonov et al.'s SAQC follows the 
prescription explained in Section |III| The unperturbed 
Hamiltonian for the AAQC is 



H Q = H B Z^ + (Hp + £ P )± (c 



(41) 



where Z^ and are projectors of a control register 
made of a qubit, and Ep is non-negative but smaller than 
1, which is the first excited energy of Hp. The ground 
and the first excited states of H are 



and 



17(0)) ® |0) = |0>' 



|»7>®|1>, 



)(n+L+l) 



(42) 



(43) 



respectively, and they are non-degenerate. Our unper- 
turbed Floquet operator is 



U = e 



-iH T 



(44) 



where T is the period of our quantum map. Since Uq 
is induced by Ho, which is composed only by few-body 
interactions, Uo can be implemented efficiently. To en- 
sure that Uo has no quasienergy between < E < Ep, 
T is assumed to be smaller than 2tt/W, where W is the 
maximum eigenenergy of Ho- From Eq. (40 1, we have 
W = E P + Wp = Ep + 2 + C(L _1 ). 

To build a quasienergy path that connects between | — ) 
and | + ) (Eqs. d42j and j43j)), we employ 



)) 



(45) 



to realize quasienergy anholonomy with iV-independent 
gap in the Floquet operator 



U s = Uoe-^W. 



(46) 



It remains to show that the kick part e ls \ v )( v \ is effi- 
ciently implementable. In the case of Grover's quantum 



search, this brings us disastrous slowdown to prepare 
\v) (Section IV A I. Contrary to this, we will show that 



the cost for the preparation of \v) is negligibly small, in 
the sense of quantum circuit complexity, to simulate the 
quantum circuit. The crucial point is to obtain \r]} (in 
| + )) from |7(0)). This is decomposed into two steps. We 



first introduce a Fourier transformation on 7i^ ck as 



^clock = l C j)( C C?)l> 
3=0 



where 



1 L 

= Y< 



-il-nkl I \L+1) 



c(0). 



(47) 



(48) 



The classical FFT algorithm ensures that the cost of 
F c i ock is O(LlnL) steps. F clock acts on |t(0)) , yielding 



lock 



l7(0)> 



L 

|a(0))®y>(Z)>. (49) 



1 

Second, we introduce a "conditional evolution operator" 



u h = [ z L +i L u L ) [ Zl^ + Il^Ul-! ; 

• ■ 1,1 } j . (50) 

which is a product of L few-body unitary operators and 
is accordingly efficiently implementable. Because its ac- 
tion on |c(Z)) extracts a product of quantum circuits as 
Uh\c(l)) = U1U1-1 ■ ■ ■ Ui\c(l)), it turns out that an effi- 
ciently implementable unitary operator U^Fdo^ gives us 
\rf) from |7(0)) 

1 L 

E/h^ciock|7(0)) = - 7 = £ |7(0> = l»?>- (51) 
VL+ 1 1=0 

Hence we have a unitary transformation G, which makes 
\v) from I — ), as follows: 



G = [Z^ + ?7 h ^ clock X (c) ) & 



(c) 



(52) 



where H^ is the Hadamard operation on the control 
register. Now it is straightforward to see that G and the 
kick operation 



-is\v){v\ _ Q e - is \-)(-\Qi 



(53) 



are also efficiently implementable. 

Once we build the Floquet operator U s , \ + ) is ob- 
tained from I — ), which is easy to prepare, through an 
adiabatic increment of s from to 2tt along the adiabatic 
passage build on Cheon's quasienergy and eigenspace an- 
holonomies of U s . Now it is straightforward to obtain 
I a), which the final state of the quantum circuit, from 

1+) im 
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To determine the running time of our AAQC, the min- 
imum quasienergy gap along the adiabatic passage needs 
to be estimated. Due to our choice of \v), the quasienergy 
spectrum of U s is simple. Two quasienergies depend 
on s. The other quasienergies are independent with s, 
because the corresponding eigenvectors have no overlap 
with \v) |16j . Among them, the relevant quasienergy is 
Ep + A, which correspond to the first excited energy of 
Hp. The distance between Ep+A and the quasienergy of 
the adiabatic passage takes minimum value A at s = 2ir, 
at the end of the passage. Hence the minimum gap is 
(D(L~ 2 ). Accordingly, the rate of adiabatic change of s 
needs to be algebraically slow in L to ignore errors due 
to nonadiabatic transitions. Thus an efficient implemen- 
tation of the simulator of the quantum circuit is shown 
to be possible. 



VI. SUMMARY AND OUTLOOK 



Adiabatic quantum computation that is originally 
composed of a slowly varying Hamiltonian is reformu- 
lated purely in terms of discrete time evolution, e.g. con- 
catenations of quantum circuits, where the concept of 
quasienergy, instead of eigenenergy, is employed to con- 
struct the adiabatic passage. Furthermore, a nontriv- 
ial family of DAQC is introduced with the help of the 
quasienergy anholonomy, which is far more easier to im- 
plement than the eigenvalue anholonomy in adiabatic 
Hamiltonian evolutions. We explain theoretical treat- 
ment of AAQC through estimations of quasienergy gaps 
for "impractical" and "realistic" cases. It turns out that 
the former case provides a key to show the equivalence 
of AAQC and the conventional quantum computation. 
The proof of the equivalence was obtained by extending 
the proof by Aharonov et al. for SAQC. Although our 
argument focus only on the essential point of Aharonov 
et al.'s proof, the rest of sophisticated points could be 
taken into account straightforwardly for AAQC. 

There remains an open question whether AAQC really 
solves classically intractable problems efficiently To clar- 
ify this, we need to examine AAQC for a hard problem by 
a classical computer. This certainly involves a simulation 
of complex quantum dynamics. In this respect, we ex- 
plain an advantage of AAQC over SAQC from a historical 
perspective obtained through the studies of classical and 
quantum chaos [28., i.e, prototypes of complex dynamics 
in classical and quantum mechanics. Although examples 
of chaos are described by differential equations [25] , the 
discovery of chaotic iterative mapping [18| [30] has galva- 
nized the studies of chaos. This approach allows us to 
facilitate both numerical and exact analysis [31 , which 
can also be applied to adiabatic quantum computation. 
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APPENDIX A: ESTIMATION OF 
QUASIENERGIES FOR GROVER PROBLEM 
WITH "FAIR" \v) 

We will examine the quasienergies of AAQC's U s 
(Eq. ([8])) for Grover's problem with the fair choice of 
\v) (Eq. (fl8|), introduced in Sec. IVB It is shown above 
that the subspace spanned by | ± ) and | /) (22 1 is the 
relevant one for AAQC. 



A perturbation expansion with a small parameter e is 
suffice to obtain the asymptotic behavior of the mini- 
mum quasienergy gap for N ^> 1. The reference of the 
expansion is a unitary operator 



W s = U e 



s\u){u\ 



(Al) 



where \u) = a\ — )+b\f) agrees with \v) as e — > 0. Because 
of the fact (u\ + ) — 0, | + ) is a trivial eigenvector of W s , 
in the sense that the corresponding eigenvalue e~ lErT is 
independent with s. Let us denote |w_(s)) and \wf(s)) 
be the two eigenvectors of W s in the remaining subspace 
spanned by | — ) and |/), where we impose |u>_(0)) = 
| — ) and \wf(0)) = \ f). The corresponding quasienergies 
W-(s) and Wf(s) are monotonically increase with s and 
exhibit quasienergy anholonomy: W_(0) = < W-(s) < 
W-{2tt) = Ep + a = W f (0) < W f (s) < W)(2tt) = 2n/T. 
At the same time, two quasienergies W-(s) and e _lBpT 
of W s crosses between < s = s c < 2tt. Hence the 
adiabatic passages along both W-(s) is inapplicable for 
the adiabatic search, because the destination of the state 
vector |/) does not provide \x) immediately. Hence we 
do need to incorporate the effect of small e to achieve the 
anholonomic adiabatic quantum search with U s 

We rearrange U s to carry out the perturbation expan- 
sion from 



where 



S s 



e +is\u){u\ e -is\v){v\ ^ 



(A2) 



(A3) 



The perturbation S s makes an exact crossing of the 
quasienergies W-(s) and e~ lEpT of W s , at s = s c , to 
an avoided crossing of the exact quasienergies E_(s) and 
E + (s) of U s . The resultant quasienergy E-(s) connects 
| — ) at s = and | + ) at s — 2ir. The narrow- 
est quasienergy gap between E-(s) and E + (s) around 
s = s c determines the speed of the anholonomic quan- 
tum search. By using the fact that both \u) and \v) 
belongs to a two-dimensional subspace of TIa ® Tic, S a 
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is expressed as S s — e lSs , where s s is the perturbation 
"Hamiltonian" 

s s x (6/sin6) _1 
= (\5)(u\ + \u)(5\ + \6)(6\)sms 

+ z(\u){S\-\S)(u\)(l-\(5\S))(l-co SS ) 

- (2\u)(u\ + \u)(6\ + \S)(u\) %((u\S)) (1 - coss) , 

(A4) 

where \S) = \v) — \u) and satisfies 

cos6 = 1 - 2 ((5\5) - \(u\S)\ 2 ) sin 2 -. (A5) 



From (5\S) = Ab 2 sin 2 § and (u\S) = -2b 2 sin 2 §, we have 

cos 6 = 1 - 8b 2 sin 2 | (l - 6 2 sin 2 sin 2 | . (A6) 

Hence, in the asymptotic regime |e| <C 1, we obtain 
e/sin0 = 1 + 0(e 2 ). Finally, we obtain the leading 
contribution of s s with respect to an expansion by e: 

s a =2e6sin| (| + )(u|e ie - ls / 2 + h.c.) +0(e 2 ). (A7) 

Thus we conclude that, if we avoid s c — 0, n, the 
minimum quasiencrgy gap between E—(s) and E + (s) is 
|2e6sin(s c /2)| = 0(e). 
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